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Abstract 

We study the effects of disorder on a system of two coupled chain of strongly 
correlated fermions (ladder system), using a renormalization group technique. 
The stability of the phases of the pure system has been investigated as a 
function of interactions both for fermions with spin and spinless fermions. 
For spinless fermions the repulsive side is strongly localized whereas the sys- 
tem with attractive interactions is stable with respect to disorder, at variance 
with the single chain case. For fermions with spins, the repulsive side is also 
localized, and in particular the d-wave superconducting phase found for the 
pure system is totally destroyed by an arbitrarily small amount of disorder. 
On the other hand the attractive side is again remarkably stable with respect 
to localization. We have also computed the charge stiffness, the localization 
length and the temperature dependence of the conductivity for the various 
phases. In the range of parameter where d-wave superconductivity would oc- 
cur for the pure system the conductivity is found to decrease monotonically 
with temperature, even at high temperature, and we discuss this surprising 
result. For a model with one site repulsion and nearest neighbor attraction, 
the most stable phase is an orbital antiferromagnet . Although this phase has 
no divergent superconducting fluctuation it can have a divergent conductivity 
at low temperature. Finally, to make comparison of our results with exper- 
imental ladder systems, we treated the interladder coupling in a mean field 
approximation. We argue based on our results that the superconductivity 
observed in some of these compounds cannot be a simple stabilization of the 
d-wave phase found for a pure single ladder. Application of our results to sys- 
tems such as quantum wires is also discussed. In particular the corrections to 
conductance in a two channel quantum wire have been obtained as a function 
of system length, temperature and interactions. 



I. INTRODUCTION 



Strongly interacting systems constitute nowadays one of the most challenging problem 
of condensed matter physics. In one dimension a fairly complete solution of the interacting 
problem can be obtained, and it is well known that one dimensional systems are one of the 
simplest realizations of non-Fermi liquids, and have generic properties known as Luttinger 



1 



liquidlhi. Prompted by a variety of experimental situations ranging from organic conductors 
to High Tc superconductors, there has been in the recent years, a growing interest in systems 
of coupled interacting electron chains. Unfortunately, despite the good understanding of 
purely one dimensional systems, the effects of interchain hopping, allowing to go from one to 
higher (two or three) dimensions are much less known. Whether non-Fermi liquid properties 
can be retained even in presence of finite hopping or not is still a highly controversial issuelli. 

Many studies have therefore focused on systems of few coupled chains! (two coupled 
chains being the so-called ladder systems), for which much more controlled analyticaluil3 
or numerical0^iil techniques can be applied allowing for a deeper understanding of their 
physical properties. For commensurate filling, i.e. one electron per site, the system becomes 
equivalent to coupled spin chains, since the charge degrees of freedom are frozen by a Mott 
transition. Important differences between ladders with and even and odd number of legs 
were expected, in a way reminiscent of the Haldane conjecture between one- dimensional 
systems with integer and half integer spins. In particular ladders with an even number of 
legs were predicted to have a spin gap. Good experimental realizations of such coupled spin 
chains like Srn-iCun+i 02nii and VOaPsOri'ii compounds have confirmed such behavior. 
Due to the presence of such a spin gap an even more spectacular effect is expected upon 
doping. At the opposite of single chains, that exhibit either a spin density wave or charge 
density wave ground state for repulsive interactions, the ladder system is believed to have 
a superconducting ground state involving pairing across the chains. That superconducting 
state has similarities with d-wave paring that has been advocated in some two dimensional 
models of strongly correlated electrons for High Tc superconductorJ^l^ such as the exis- 
tence of a spin gap and a sign change of the superconducting order parameter when one 
moves on the "Fermi Surface". In the strong coupling limit i.e. the t-J model, the d-wave 
phase can also be viewed as a RVB stateil. 

However all the studies of ladder systems have been, up to now, restricted to pure sys- 
tems. Unfortunately (or maybe fortunately) it is well known, that for one dimensional 
systems disorder has extremely strong effects. For a non-interacting system, it is well known 
that all states get localized in the presence of an infinitesimal random potentiaSi. In- 
teractions can modify this picture, but for one chain system delocalization occurs only for 
strongly attractive interactions. In particular even normal s-wave superconducting phases 
are destroyed by non magnetic impurities except for exceedingly attractive interactions (see 
e.g. Ref. |36| and references therein), and no Anderson's theorem exists even for weakly 
coupled one-dimensional systems^. 

In order to compare the theoretical predictions of d-wave superconductivity in doped 
ladder systems with experiments, it is therefore of prime importance to understand the 
effects of disorder on the phase diagram of the pure ladder system. One of the important 
question is of course the stability of the newly found d-wave superconducting phase since 
there are no obvious reasons why it would survive the introduction of a small amount of 
disorder. Such a study is also relevant to the physics of quantum wires with few channelS^. 
In quantum wire systems, the situation is however complicated by the occurrence of long 
range Coulomb forces0 that can induce a one dimensional analog of the Wigner Crystal@'0, 
which can drastically modify the response of the system to disorder^. However, the presence 
of charges in the grids of the quantum wire systems can be cleverly used to screen completely 
the long range interactions, and have an experimental realization of a Luttinger liquidcil. By 
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changing the gate voltage it is possible to have more than one band at the Fermi level, in 
a controlled way. The quantum wire is thus a possible realization of a two (or more) leg 
ladder. Interband tunneling plays the role of interchain hopping. They provide ideal systems 
in which to check for the effect of disordeiMEll. 

Besides the exciting possibility to test the ability of Luttinger liquid models to describe 
accurately the now available quasi one dimensional experimental systems, investigation of 
disorder effects in ladders presents in its own right a great theoretical interest. Indeed, the 
two chain problem is the simplest one to study the effects of interchain hopping onto the 
Anderson localization in presence of interactions, giving some clues onto this difficult topic 
in more than one dimension. In particular, one would be interested in obtaining boundaries 
between localized and delocalized phases and the dependence on localization lengths on 
disorder. Another question of particular interest is the effect of interactions on physical 
quantities controlled by disorder such as the conductivity for a macroscopic system, or for 
a mesoscopic one the persistent currents. In particular, for a one chain system, it was 
shown that for a system with spin degrees of freedom persistent current were enhancecH by 
repulsive interactions, at variance on what happened for a spinless system. It is therefore 
important to check whether this striking results still holds in a more two dimensional system. 

In this paper, we consider the effects of a weak random potential scattering on systems of 
coupled fermionic chains both with spin and spinless using bosonization and RG techniques. 
A short account of some of the results of this paper were presented in Ref. Besides giving 
the phase boundaries, the RG also provides us with expressions of the localization lengths, 
temperature dependence of conductivity, and dependence of persistent currents with system 
size. The plan of the paper is as follows : 

In section |I|, we discuss the spinless fermions 2 legs ladder problem . We first recall the 
phase diagram of the pure systemlll, then consider the effects of disorder. This allows for a 
detailed comparison of the transport properties of the 2 chain system with the ones of the 
one chain spinless fermion system and the ones of one chain of fermion with spin. We show 
that contrarily to naive expectations, the ladder spinless fermions system is very different 
from the one chain system with spin, and that the effect of interactions on persistent currents 
is even more violent on a two chain spinless fermions system than it was for a one chain 
system. 

In section |T|, we discuss the technically more involved case of fermions with spins. 
Following the same methodology, we first recall the phase diagram of the non-random 2 
chain system and then consider the effects of a weak random potential on the phase diagram. 
As for the spinless problem, we give a detailed discussions of the transport properties in 
the disordered phases. We compare these results with the ones already known for one 
chain, and show that the reduction of persistent current by attractive interactions that 
is observed in one chain of fermions with spin should be almost absent in the two chain 
system. The d-wave superconductivity of the fermionic 2 chain system is a feature that is 
not preserved in the presence of a very small amount of disorder. On the other hand for 
some values of the parameter, and orbital antiferromagnetic phase exists. This phase has an 
infinite conductivity even in the presence of disorder although it has exponentially decaying 
superconducting correlations. 

In order to compare our results with experiments on doped ladder systems, it is necessary 
to treat interchain coupling which stabilizes superconductivity at finite temperature in real 
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systems, and may also reduce the sensitivity of the system to disorder. Thus in section 
IV we examine the mean field theory for the d-wave superconductor in an array of coupled 
disordered ladders. We give a criterion for persistence of superconductivity in the presence 
of disorder and show that d-wave superconductivity remains unstable except in very pure 
samples or in the presence of a very strong Josephson coupling between ladders. 

In section 0, we summarize the implications for experimental systems such as doped 
SrCuO chains and quantum wire with two channels. We claim that in recently synthesized 
doped ladder systems, the physics of the superconducting phase is more likely to be of two 
dimensional origin rather than just a stabilization of a ladder d-wave superconductivity. 
Conclusions can be found in section |V^. Finally most of the technicalities can be found in 
the appendixes. 



II. SPINLESS FERMIONS 



A. Pure system 



Let us consider first two chains of spinless fermions coupled by an interchain hopping t± . 
For simplicity we first consider only nearest neighbor interactions, and interchain interac- 
tions. The effect of more complicated interactions will be detailed below. The Hamiltonian 
for the pure system reads 



H = -tY^ 4pQ+i,p + h.c. + VY^ rii^pUi+i^p 



+t± J2 4,iCi -1 + h.c. + UJ2 n'i,ini_i 



(1) 



where p = —1,1 is the chain index and i is the site index. To treat the interactions, it 
is convenient to rewrite the Hamiltonian in term of boson operatorsi'Sil. To do so, we 
linearize the fermions dispersion relation around kp, introduce right (R) and left movers (L) 
for each chain, and take the continuum limit Cn^r,p y/(^'4'r,p{na) with r = L, R, p = ±1 the 
chain index and a the lattice spacing. We use the bonding tpo = and anti-bonding 

V^TT = bands base and introduce the densities Pr,o,iT{x) =■ ipl^g^^{x)ipr,o,iT{x) :. We then 

define the canonically conjugate fields 0pj| and Hp y via 



dx(l>p,\\ = ^'^^PL,o + PR,o ± pL,n ± PR,7t) 
np,|| = -^iPR,o ± PR,7T - pL,o =F pL,n) 



(2) 
(3) 



and the field 6n 



J^^Ilpn{x')dx' . More details on the bosonization technique can be 



found in Appendix 0. In term of these fields the Hamiltonian becomes 
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H = Hp + H\i 



dx 
2^ 



(4) 



m= I- 

27r 



«,|7r||(7rn||)2 + ^(9.0ll)2 
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For the microscopic Hamiltonian (|l|), one finds 
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(5) 



with vf = 2tasm{2kFa). Therefore for the pure t-V model, one has Kp < 1 (resp. Kp > 1) 
and Qf < (resp. Qf > 0) for repulsive (resp. attractive) interactions and K\\ = 1 for all 
t, V. In fact (^ describes the most general two chain spinless system. More complicated (i.e. 
longer range and interchain interactions) lead only to a change in the parameters K, u and 
g. By adding interchain interactions such as U in formula (|^) or longer range interactions 
one can in particular access the other regimes Kp > 1 and Qf < or Kp < 1 and Qf > 0. The 
physics of the system is readily seen on (^). The t± term suppresses cos(v^0||). Depending 
on the value of K^\, the can either remain massless or develop a gap. We concentrate here 
on the case where develops a gap and acquires a non-zero expectation value determined 
by minimizingthe ground state energy (see appendix ^). This situation always occur for 
the t-V modelu. By mapping (^) on a problem of one chain of fermions with spin and spin- 
anisotropic interactions in a magnetic field^, one can obtain the complete phase diagram for 
the pure caselli. Since, due to the one dimensional nature of the problem, no true ordered 
state exists, one has to find the most divergent instability. As for one chain, two main type of 
instabilities are possible: particle- hole (density, current etc.) instabilities or particle-particle 
(i.e. superconducting) ones. The operators with the most divergent susceptibilities are in 
boson form 

OcDW^ = <i(x)^L,i(x) - ^lj^_i^i,-i(x) ~ e*^'^''cos(v^e||), 
Osc^ = i^L,o{x)i^R,^ + i^L,n^R,o ~ e*"^^" sin(V2^||) 
OoAF = t{^l,{x)^L,^,{x) - tpl_,{x)^L,,{x)) ~ e^^^" sin(v^^|,). 



~ e*^^" cos(V2^||) 



They describe respectively out of phase charge density waves, an orbital antiferromagnetic 
phase and chain symmetric "s" and antisymmetric "d" type superconductivity. The out of 
phase charge density has a 2kF modulation of the density along the chain and a change of 
sign across the chains. In the orbital antiferromagnet currents go from one chain to the 
other with wavevector 2kF, giving currents circulating around plaquettes of length n/kF- 
The superconducting phases are the standard ones, given on the original model by 
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OsCd{n) = Cn,lCn,2 



(6) 



The most stable phase depends on the parameters K and g. The various cases are given in 
Table |, and the phase diagram shown in figure |I]. In Ref. |T3| the bosonized forms of Osc" 



and Osc<i are exchanged due to the neglect of anticommuting operators (see appendix 0), 
so that the two superconducting phases have been erroneously exchanged. 



B. Effects of disorder 



Now, we consider the effect of the disorder on (|T]-^. We introduce a random on-site 
potential e^^p uncorrelated from site to site and from chain to chain: 



-^random ^ ^ ^i,p^i,p^i,p i.'^^ 



P=±l 



with ei^pCj^p' = DSijSp^p'. In the continuum limit and using the bonding antibonding basis 
the disorder becomes 



-^random / dx 



es{x){i'lix)'ipo{x) + ipl{x)ip^{x)) + ea{x){'4)l{x)i)^{x) + V^;^ (x)^o(a;)) (8) 



with es,a = (ei ± e-i)/2 and eQ,(x)e/3(x') = ^S{x — x')Sa,i3- Using the expression of fermion 
operators defined in Appendix ^ and (0) one obtains for the disorder term 



i^random = [dx Vs{x)—d,M^) + ^iMg^v^^P COs(v^0,|) + ^Mg-V^-^P COs(v^0|| 

+ [dx cos(v^0i|) cos(v^^^ii) + ^e'^'^^ cos(v^0||) + ^Mg-^v^^p cos(v^0,| 



(9) 



where the disorder has been split in a g ~ component {r]s,a) and a g ~ 2kF one (^s,a)- As 
for one chain the rj and ^ are uncorrelated and 



Vs,a 




= Ds,aaS{x - x' 




{x)^s,a{^') 


= 



(10) 

(11) 

6,a(x)e,a(a;') = Ds,aa5ix - x') (12) 

The g ~ (forward scattering) part of the disorder does not affect the conductivity 
and cannot lead to localization^, but could in principle modify the phase diagram and in 
particular destroy the gaps of the pure phase. As for one chain, one can eliminate the 'r]sdx4>p 

by a transformation (pp (pp + J dxr]s{x). The only effect of this term is therefore to 
give an additional exponential decay in the density-density correlation functions. 

Due to the presence of a gap in 6*11 (see table D, the T]a{x) cos(a/20||) cos(-\/26'||) term is 
always suppressed at lowest order. It could however generate relevant terms at higher order. 
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However higher order terms are either identical to backscattering terms aheady present in 
the Hamiltonian, or adds random contributions to cos a/80 y and gf cos y/SO^i. At small 
disorder these contributions are negligeable, and one can completely disregard the forward 
scattering. We can therefore keep only for the coupling to disorder Hs + Ha 

Hs= f —Ux)e'^^^ cos(V2</.||) + h.c. (13) 

Ha= f —U^)e'^^'' cos(v^^^ii) + h.c. (14) 

In Hs the symmetric part of the disorder couples to the in-phase charge density wave order 
parameter Ocdw^ = cos(\/20||), whereas the antisymmetric part involves Ocdw^- 

Due to the gap in 9\\, (f)\\ has huge quantum fluctuations, and consequently the symmetric 
part of the disorder Dg is always less relevant than the antisymmetric one Da- We can 
therefore focus on the latter and forget about the former. The effect of (|T^) again depends 
on the values of gf and K. 



1- 9f<0 

For gj < (i.e. ^ > for the t-V model) we can replace cos{^/29^\) by its (non-zero) 
mean value and the coupling to disorder ( p!^ reduces to C / dxC,a{x)e^^'^''^^'' + h.c, where C 
is a constant. The effect of such a term can be determined, as for a single chainS, by using a 
renormalization group (RG) procedure. Upon varying a cutoff a, similar to a lattice spacing 
in the original lattice problem, one find the following renormalization for the disorder 

^ = B,(3 - A'p) (16) 

where / = ln(a) and C2 a constant. (^) implies a localization-delocalization transition at 
Kp = 3. For Kp > 3 the disorder is irrelevant and the corresponding phase in the pure system 
is stable. For Kp < 3 disorder grows. Although the system flows to a strong coupling fixed 
point, it is naturalil to interpret this phase as localized by disorder, since the disorder will 
pin the massless field (j)p. As a consequence, the d-wave superconducting phase is unstable 
in the presence of disorder except for huge attractive interactions. In the case of the t-V 
model at \^ > 0, we have Kp < 1 and therefore the CDW^ is always pinned by the disorder. 

Similarly to the one chain problem the localization length can be computed using the 
RG. For very weak disorder and far from the transition one can neglect the renormalization 
of the exponent Kp induced by Da- Using that approximation, we obtain: 

Da{l) = e(3-^''')'D„(0) (17) 

For Da{l) ^ v\la that scheme breaks down and we have a strongly disordered system. For 
such a system the localization length, i.e. the scale of variation of the phase 0p is of the order 



of the (renormalized) lattice spacing a* . This occurs for e'* ~ ( ^ ^^^^^ ) . Therefore 
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/ ,,2 \ 3-Kp 

Let us recall that for a non-interacting system, the localization length is of the order of the 

mean free path i. e. Lioc. ~ -pa- 
using the renormalization equation it is also possible^ to obtain the temperature de- 
pendence of the conductivity for temperatures above the the pinning temperature U/L2 ch.- 
Below the pinning temperature, the conductivity is expected to decrease as exp — (Tpi^jT)^, 
by analogy with non-interacting electrons. A derivation of the temperature dependence of 
conductivity has been given in Ref. Another method to derive the temperature (or 
frequency) dependence of the conductivity is given in appendix |B[ If one neglects the renor- 
malization of the exponents the conductivity behaves as 

a{T) ~ T^-^" (19) 

Therefore, for Kp < 2, the conductivity decreases, and there is no remnant of any su- 
perconducting behavior effect well above the temperature at which the system is effectively 
pinned Tpin ~ u/L2ch.- Thus the existence of d-wave superconductivity in the pure system 
affects the transport properties of the disordered system only for quite large attraction. 
Analogous effects will occur for fermions with spins as will be discussed in section |IT1[ 



2. gf>0 



For gj > (i.e. attractive interactions for a t-V model), (6*11) = and in a first 
approximation the coupling (n) vanishes. Obviously, this approximation is too crude and 
one must integrate the fluctuations of 6\\ around its mean value to get the effective coupling. 
This is done in appendix and gives the following effective action for (pp-. 



dxdr 



(V0p) 



+ h.c. 



(20) 



with i{x)i*{x') = DS{x - x') and D ~ D^. 

The renormalization of the disorder is given by an equation similar to (P^: 



dD 
111 



{?,~AKp)D{l) 



(21) 



The disorder is now relevant only for Kp < 3/4, leading to three different phases for gf > 0: 
a random orbital antiferromagnet for Kp < 3/4, an ordered orbital antiferromagnet for 
3/4 < Kp < 1 and a s-wave superconducting phase for Kp > 1. For the t-V model, 



Kp > 1, and the 



s"-wave superconducting phase is therefore stable with respect to weak 
disorder, at variance to the single chain problem. For the latter the delocalization only 
occured for extremely attractive interactions i.e. Kp > 3/2. For the two chains problem the 
localization-delocalization transition arises in the immediate vicinity of the non-interacting 
point. Contrarily to the case of repulsive interactions, interchain hopping now strongly 
reduces the localization effects. 
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The localization length in the random orbital antiferromagnet, is now given by 



2 



The conductivity behaves both in the OAF and the s-wave phase as 

a{T) ~ T^-^^" (23) 

diverges as T since the ground state is superconducting. It is to be noted that although 
the OAF has no superconducting order parameter, its conductivity can also be divergent 
for Kp > 3/4 even in the presence of disorder. An expanded discussion of Orbital Antifer- 
romagnet phases can be found in [111 C 1| and appendix 0. 



The resulting phase diagram is summarized on figure ^ together with the single chain 
phase diagram. 



C. Physical consequences 

The ladder system shows drastically different sensitivity to disorder depending on the 
sign of 5^/: at (7/ < localization effects are much stronger than at (?/ > 0. This is obvious 
both on the phase diagram shown on figure ^ and on the expression ([T8|) and (|2^) for the 
localization length. For the case of a pure t-V model, (?/ > Kp > 1 when V < (attractive 
interactions) and as can be seen from figure the system is delocalized. Although our 
calculation do not allow us to come arbitrarily close to the V = point for finite disorder, 
since the disorder has to be smaller than the gaps of the pure system, we see that if we have 
a very small disorder, the insulator superconductor transition does occur in the vicinity of 
the non-interacting point. This is remarkable and in marked contrast with the single chain 
system where the delocalization transition occurs for K = 3/2 i.e. very strongly attractive 
interactions even for arbitrarily weak disorder. One could naively think that this effect is 
simply a manifestation of the delocalization effect seen for non-interacting electrons when one 
increases the number of channel (or the number of chains). The mechanism is more subtle 
however, and in in fact controlled by the interactions. Contrarily to the noninteracting 
case where the localization length is simply proportional to the number of chains, we have 
here a complete delocalization of the attractive region, and the localization length becomes 
infinite. 

For the repulsive case V > (i.e. (?/ < Kp < 1) the opposite effect occurs and ladder 
system is more localized than the corresponding one chain system. Indeed for one chain the 
localization length is given bySS 



Li ch. / 



a \ Da 



(24) 



and is therefore longer than the one of the ladder system shown in ([T^). For very large 
repulsion {K 0) these two le ng th give back the standard Fukuyama-Lee pinning length 
of classical charge density waves^. For finite repulsion the localization length of the ladder 
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system is much shorter than the one of the corresponding one- dimensional system with the 
same K. Close to the non interacting point ~ 1, the localization length of open chain is 

just the mean free path Li ch. ^ Vp/D whereas the ladder one is L2 ch. ~ « 

This peculiar behavior of the spinless ladder system is due to the gaping of some charge 
modes, that is different depending on whether the interaction is attractive or repulsive. For 
the repulsive side 2kF charge fluctuations are still there and the gap just reduces some of the 
quantum fluctuation and hence reinforce the effects of disorder, whereas for the attractive 
side the gap kills the dominant charge fluctuation coupled to disorder and helps to delocalize. 
The sensitivity to disorder is therefore not directly related to the presence or absence of 
superconducting fluctuations in the pure system, but more on how the density fluctuations 
behave. The smoother are the density fluctuations, the less localized the system is. These 
effects will be even more transparent for the system with spins as will be examined in details 
in section As a consequence the transport properties cannot simply be guessed by looking 
at the phase diagram of the pure system. They even can be opposite to what our intuition 
based on higher dimensional system could suggest: the more "superconducting" the system 
is the better is the transport (see e.g. section |111 C| ). 




D. Persistent currents in the ladder system 

In addition to the temperature dependence of the conductivity, one can compute the 
charge stiffness of the system!!^ D, which measures the strength of the Drude peak in a 
macroscopic system a{uj) = V6{uj) + a^eg- The stiffness V can be related to the change of 
the energy of the ground state of the system in presence of an external flux by 



L £Eq 
V = 

2 dd)'^ 



(25) 



0=0 



Eq being the ground state energy of a ring in a fleld. denotes the boundary angle = 
27r///o where / is the flux threading the ring and /o = hc/e is the flux quantum. This 
quantity is directly related to the persistent currents for a mesoscopic systemS^ For a 
mesoscopic system, the persistent current measures the response to a flnite flux by 



(26) 

Therefore the stiffness T) provides a measure of the persistent currents for small (or close to 
a multiple of 2tt) flux since J = 2T>(f). Although the complete calculation of the persistent 
currents at flnite flux is also possible for a one dimensional interacting system, the calculation 
is more complicated in the presence of disorder, and the stiffness carries enough information 
for our present purposes. 

The effects of interactions on persistent currents is an extremely difficult question to an- 
swer in two or three dimensions. Perturbative calculations suggests that interactions could 
enhance persistent currentsEHH^. For a single spinless chain the persistent currents were 
found to decrease with more repulsive interactional^. This effect can naturally be ex- 
plained using a renormalization group technique, and it was shown that such behavior is 
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peculiar to the spinless problem and that for a single chain of electrons with spins persistent 
currents should be enhanced by repulsive interactions@@. For the ladder system it is there- 
fore very interesting to see if the same effects occur and in particular to check again for the 
differences between the spinless system and the system with spins. In particular one could 
imagine that the chain index acts in a similar way than a spin index for a single chain. As 
we will see this idea is far too naive. We examine the spinless system in this chapter and 
the system with spins will be investigated in chapter |ITT|. 

For the ladder system, the conductivity stiffnessSis obtained using ( |A13| ) as "D = 2upKp. 
The factor of two compared to the single chain expression ( [A13| ) is due to the fact that there 
are twice as much degrees of freedom in the 2 chain system. In the following, we consider a 
finite system, the size L of which is smaller than the localization length. 

From the renormalization group equation for Wp, one can obtain@ the renormaliza- 

tion group equation for V 

The conductivity stiffness of a disordered system of size L, 'D[L) is then obtained by stopping 
the RG equation at a{l) = L and taking T>{L) = T>{1). In the case Qf < 0, we have seen 
that D{1) = D(0)e^^^^''^'' , at least when a{l) <^ L2 ch.- Putting that approximation for D{1) 
in (p7|) gives us: 



V{L) = V{0)-CD{0) 



L 



3-K„ 



.«(o). 



(28) 



Using the expression for L2ch., (0) simplifies for length smaller than the localization length 
into into 

1 -1 



P,,<o(L) = P(0) - C 
P,^,>o(L) = V{0) - C 



^2ch.{gf<0), 



'2ch.igf>0), 



3-4Kp 



(29) 
(30) 



Thus for gj > the reduction of the stiffness is less important than for Qf < 0. 

Therefore, the length dependence of the conductivity stiffness (and the persistent cur- 
rents) is extremely sensitive to the attractive or repulsive character of the interactions for 
the t-V model or any model with intrachain-only interactions. By comparison with the one 
chain case@@, we see that the effects of the interactions on the conductivity stiffness are 
qualitatively the same (i. e. repulsive interactions help in reducing the conductivity stiff- 
ness, while attractive interactions reduce the decrease of conductivity stiffness by disorder) 
but they are much stronger for two chains than for one chain. In fact, for a t-V model, 
the reduction of conductivity stiffness would be finite for attractive interactions, even in an 
infinite system since then the disorder is completely irrelevant. 

It is noteworthy that the chain index does not act in a similar way as a spin degree of 
freedom, for which there would be an increase of the persistent currents showing again the 
important difference between a system with and without spin@. The physical reasons for 
this difference are examined in more details in the next section. 
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E. Spinless ladder vs. one chain with spin 



Naively, one could think that going from one chain to two chain amounts to having 
one internal degree of freedom that is equivalent to spin, and thus that the results for the 
system with spin will apply straightforwardly to the ladder system. However, from what 
we have seen precedingly, this is definitely not the case. In fact, we have properties for the 
spinless ladder that are just the contrary of the ones of the fermions with spin. Attractive 
interactions delocalize in the spinless fermions case, whereas they increase localization in the 
case of fermions with spin. Persistent currents are enhanced for more attractive interactions 
in the spinless ladder whereas repulsive interactions would enhance the persistent currents^ 
in a spin system. The reason for that is that the spinless ladder has no SU(2) symmetry 
(except for V = 0) contrarily to one chain with spin. The minimum of the ground state 
energy of the spinless ladder corresponds to states that break the SU(2) symmetry because 
t± plays the role of a magnetic fieldll3c2l. Thus such phases cannot be obtained in an isotropic 
system of fermions with spin. 

For attractive interactions, the only way for the symmetric fermions with spin system 
to preserve SU(2) symmetry is to form singlet phases such as 2k p charge density waves or 
singlet superconducting state. Coupling the charge density wave fluctuations with a random 
potential implies strong localization effects. On the other hand, the spinless ladder simply 
form pairs along the chains and can avoid to form 2kp fluctuations. Translated in the spin 
language, such phase would be an anisotropic triplet superconductor with a spin gap, and 
would be forbidden by symmetry. In the same way, for repulsive interactions, preserving 
SU(2) symmetry prevents the formation of a gap, whereas a gap formation is possible for 
the spinless ladder giving an out of phase charge density wave. In the spin language, this 
corresponds to an anisotropic SDW. 

Adding random potentials to the spinless ladder results in a rather artificial model of 
fermions in a random potential and a random field parallel to the z axis. Because of the 
anisotropy, the system is more sensitive to the random field parallel to the z axis than to the 
random potential. Thus, for repulsive interactions, the anisotropic system has a very strong 
coupling to disorder, whereas for repulsive interactions, it is only weakly coupled. On the 
other hand, the isotropic system is only feeling a random potential. When interactions are 
attractive, there is a spin gap and CDW fluctuations that can couple to disorder, making 
the system more localized. When interactions are repulsive, on the other hand, there is no 
spin gap thus reducing the coupling of the CDW fluctuations with disorder. 

We conclude that for interacting systems, contrarily to their non-interacting counter- 
parts, not only the number of available internal degrees of freedom but also the internal 
symmetries determine the response to random perturbations. Loosing some symmetries 
allows for a larger variety of ground states, and thus to very different responses to weak 
perturbations. 

III. FERMIONS WITH SPIN 
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A. Pure system 



The piire case has been analyzed in great details both analyticallyiHli and 
numerically0lii0iil. A very interesting feature of that model is the existence of a "d-wave" 
superconducting phase for purely repulsive interactions and the existence of a spin gap. The 
Hamiltonian is in the extended Hubbard case: 

^ ^ 5Z ^i+l,<J,p^h^<P "I" H. C. t± Cj^^^pCi^fj-p 

i,cr,p i,c,P 

+U ^ Tij.t.p^-ia.P + ^ H ni^pni+i,p (31) 

i,p i,p 

with p = ±1 is the chain index and a =|, | labels the spin. In order to treat this Hamiltonian 
using bosonization one has to separate the bonding o and antibonding vr bands as was done 
for spinless fermions. Then, within each band, one can apply the standard bosonization 
formulas for fermions with spins. As a consequence, the system is described by 4 fields 
0p) 0^0°, 0° instead of 2 in the spinning case. For the pure case we follow closely the 
derivation of Ref. |l^. It is convenient in the following to replace the fields (f)°''^ {v = p, cr) 
by linear combinations: 0y± = -^{(pij,o ± <Pu,Tr) The low energy physics depends on the signs 
of two constants gi,g2- Physically, g2 represents the forward scattering interaction, while gi 
represents the backward scattering interactions. The Hamiltonian consists of a free part : 

H= UurK^riTTll^rr + ^{d^.r?] (32) 



and two sine-Gordon like part, one associated with interband processes induced by intrachain 
forward scattering : 

g2 f 

Hint,2 = —7 — / dx COS 26 p_ (cos 20o-_ + COS 26'o— ) (33) 



2(7r«)2 . 

The other associated with the intrachain backward scattering : 
2(7* r 

Hint,i = 7 / dx [cos 2(/)o.+ (cos 26'p_ + cos2(f)a-- + cos26'o— ) — cos20p_ cos20o-_] (34) 

[2na) J 

In all cases, only one of the four bosonic fields (0p+) is gaplesJii and all physical quantities 
depend on a parameter Kp+ of the symmetric charge mode, analogous to the Kp of the 
spinless problem. In terms of gi,g2,Kp^ is given by : 

That expression is valid for the generic g-ological model. For the extended Hubbard model, 
we can go further as gi,g2 can be expressed in terms of U,V, kp as : 

gi = Ua + 2Vacos(2kFa) 
- 2c/2 = -{Ua + 2Va{2 - cos{2kFa)) (36) 
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where a is the lattice spacing. The mean values of the three other fields are determined by 
minimizing the energy of the ground state. Depending on the interactions one can distinguish 
four sectors that are summarized in table |1| 

As for the spinless case one has to consider the various operators with divergent suscep- 
tibilities 

OcDW-{n) = ^Pcia^pCn,a,p (37) 
p,a 

OoAF{n) = Y.P(^i,a,pCn,a,-p (38) 
OscijL) = XI Cn,<7,pC„ _(^,p (39) 

p 

OsC^i^) ='^Cn,a,pCn-a-p (40) 
P 

When taking the continuum limit these expressions become 

OcDW- = J2i'^Lla^Rl'7 - ^LlaV'R-la) (41) 

OoAF = ^ ^KV^il^rV^fl-l'T - i^{-lai^Rla) (42) 
cr 

OsC = J2 ii^LOai'm-a + IpL^alpRTT-a) (43) 
cr 

OsC-i = ii'LOa'il^RO-a " ^Lnai^RTr-a) (44) 
(J 

where for the SC operators, one has to retain the g ~ component, whilst for the OAF 
and CDW" the q ~ 2kF component gives the dominant contribution. To get the correct 
bosonized expression one has to pay extra care to the anticommuting U operators^ and one 
obtains 

(45) 
(46) 
(47) 
(48) 

From the bosonized form of these operators (simplified by the fact that {Op-) = 0) everywhere 
and the expressions given in table || one can deduce that sector I is a SC^ phase, sector II 
an OAF phase, sector III a SC^ phase and sector IV a CDW^ phase. The phase diagram 
of the pure system is summarized in figure ^. Note that for the pure Hubbard model, which 
corresponds to 1/ = in (p6D, one can only have the SC^ phase (for U > 0) or the SC 
phase (for U < 0). The other phases could be obtained for a more general model such as 
the extended Hubbard model. We will come back to that point later. 



OcDW^ 


2 

vra 


6"^"+ COS i 




sin 9a. 
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B. Effects of disorder 



Let us now add a weak random on-site potential : 

-^^random potential — ^ ] ^i,p'^i,a,p ('^^) 



l,(T,p 



with rii^p — 4-|^Cj^-f + c\^^Ci^i and Ci^^e^y = DSijSp^pi. We go through the same steps as in 
the spinless fcrmions section. We got to the continuum Umit, introduce the bonding and 
antibonding band, and bosonize the resulting coupling to disorder. Let us first consider the 
q ^ part of the coupling to disorder. For the symmetric part of the disorder this coupling 
is of the form : 

Hs,q^o = j 'ns{x)dx(j)p{x)dx (50) 

It is clear that this part of the disorder can be eliminated by the transformation (t)p{x) 
(ppix) + ^^^^rjs{x')dx' For the g ~ part of the antisymmetric random potential, we 
obtain : 

^a,q~0 = / dxr]a{x) [i^R,Q,a'^R,'n,a + '^L,0,a'^L,^,<T + H. c] (51) 

The bosonized form of that operator is the following: 

i^a,,~o= /dx^[e^('^''-+^''-)cos(0,_ + ^._)+e^(-'^''-+^''-)cos((/),_-^,_)+ H. c.l (52) 

From that equation, we see that the g ~ part of the antisymmetric disorder is not coupled 
to the gapless charge symmetric mode. Moreover, it always contain one term that has 
exponentially decaying correlations. Therefore, it cannot break any gap by an efi^ect d la 
Imry Ma and cannot generate any relevant term by a massive mode integration. It will thus 
be possible to drop it safely in the following. Then, we have to consider the 2kF part of the 
disorder. We have for the 2kF coupling to disorder two terms : 

^0, = j ^a{.x)OcDW^{.x) + Cai.^)OcDW-i.^)dx (53) 

Hs = J ^x)OcDwo{x) + Cs{x)Oljjwo{x)dx (54) 



Where ^nix)^n'ix')* = Dn6n,n'Six — x'){n,n' = a,s), the being random Gaussian dis- 
tributed potentials. The operators Ocdw° represents the in-phase charge density wave, and 
OcDW^ the out of phase one. 

As before we assume that the disorder weak enough not to destroy the gaps in the 
system. We have already argued that the g ~ is irrelevant to our problem. Concerning 
the 2kp part, we only retain the masslcss mode. The situation is quite similar to the one 
of a XXZ spin chain in a random magnetic field. A XXZ spin chain is a Hubbard chain 
at half filling, and thus has a charge gap. The random magnetic field couples to the spin 
density that contains (frozen) charge degrees of freedom. However, the random magnetic 
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field only affects the spin degrees of freedom and does not break the charge gap. By analogy, 
we expect that even when the random potential gets relevant it will not break the spin gap 
or the gap in the antisymmetric charge mode. Since the gaps are stable, we can obtain 
simplified forms for the couplings by replacing the fields by their mean values as we did in 
the spinless fermions problem. 



1. SC^ sector 



We want to analyze the effect of the weak random potential introduced through (|53| - p4D . 
Making use of the full expressions of Ocdvk°.'^ and replacing the gapped fields by their mean 
values (see sector I of table HD we obtain the following simplified forms : 

OcDvyo -6*^"+ sin(0^_) (55) 
OcDW- ~ 6*^"+ sin(e,_) cos 0,+ (56) 

These two operators have exponentially decaying correlation functions and no direct coupling 
with disorder would exist if one just took into account the mean values of the fields (pp- 
and Oa~- As in the spinless case one should integrate over fiuctuations to get the effective 
coupling 



cdisorder 



J ^esXx)e''^''''+'''''^^dxdT + H. c. (57) 



(^) can be viewed as the coupling of the fermions with the 2{kpo±kpT^) Fourier component of 
the disordered potential i. e. to a Akp charge density wave. The origin for such a 4/ci? charge 
density wave can be understood in simple terms: at half filling, the strong on site repulsion 
puts one fermion per site, meaning that there are no 2kF CDW fiuctuations. However, the 
fermion density is maximum on the lattice site and minimums in between giving the ^kp 
charge density wave fiuctuations. In addition due to the spin gap occuring in a ladder with 
an even number of legs there are no 2kp fiuctuations in the spin density as well. As we move 
away from half filling, the spin gap will survive as well as the absence of 2k p fiuctuations. 
Therefore, a random potential can only couple to the Akp component of the fermion density 
even away from half filling. This is to be contrasted to the case of a single chain where the 
dominant coupling occurs through the 2k p charge fiuctuation. One thus expects the disorder 
effects to be weaker in the ladder system. One can also recover directly the Akp CDW by 
looking at higher Fourier components of the density in the bosonization formulas. The 
physics of the metal insulator transition can be here interpreted as the pinning-depinning 
transition of this 4:kF charge density wave. 

Due to the presence of the gaps, the problem has in fact been formally reduced to a 
problem of one chain of spinless fermions with disorder . Using the results from the one 
chain problem we find that the localization-delocalization occurs at Kpj^ = 3/2. Since purely 
repulsive interaction imply Kpj^ < 1 the ci-wave phase is therefore unstable to arbitrarily 
weak disorder. The symmetric (|5^ ) and the antisymmetric (^3j) part of the disorder con- 
tribute equally to destroy the d-wave superconductivity, in contrast with the spinless case 
where the antisymmetric part was the most relevant. The localization length of the 2 chain 
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system with spin and purely repulsive interactions can be obtained by a similar method than 
for the spinless case and is 



/ „2 \ 2/(3-2Kp+) 
>2 ch. / Vp \ 



(58) 



and therefore longer than the corresponding one for one chain with repulsive interaction^^*"^^^' 

Li ch. 



a 




(59) 



As for the spinless case (|58|) is applicable if one is far enough from the noninteracting point 
so that disorder does not destroy the gaps created by the interactions. In that case one 
sees from ( p8|) that there is a considerable delocalization in the ladder. Indeed for weakly 
repulsive interactions Kpj^ ~ 1, the localization length becomes much longer than the mean 
free path /, since L2 ch. ~ a{l/af', instead of L ~ / for a single chain. However the more 
repulsive the interactions become, the more the system localizes (one recovers L2 ch ~ ^ for 
K = l/2). 

The temperature dependence of the conductivity can be obtained above the pinning 
temperature Tpin. = (see appendix g). One gets 

a{T) oc T^-^^''+ (60) 

For Kpj^ < 1, the conductivity decreases as T — even for temperatures much higher 
than Tpin . There is no remnant of the "superconducting" behavior of the pure system in 
the whole SC^ sector {Kp^ < 1). 



2. SC sector 

For sector III, the O operators take a different simplified form, due to the different gaps 
in the system 

OcDW ~ 6*'^"+ cos(0p^) (61) 
OcDW^ ^ e"^^+ sm{e^-) (62) 

By substituting in ([53| ) and (|^) and integrating over fluctuations we end with an action 
of the form (57). This time, i^p+ > 1, so the localization-delocalization transition can be 
reached at Kp^ = 3/2. This transition arises for much weaker attraction than in the one 
dimensional caseil where Kp = 3. This critical value of K can be realized for a simple Hub- 
bard model (the maximum K for the Hubbard model is i^' = 2l^3c3) whereas the one chain 
Hubbard model is always localized even for very negative f/El. In addition the localization 
length is increased 

2 

L2 ch. _(V^\ ^gg^ 



a \ Da 
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whereas in the one chain case 

= (^)3^ (64) 
a Da 

Note that here the locahzation length has the same dependence in disorder on the attractive 
(|63|) and the repulsive (|5^) side, whereas for a single chain the localization length is reduced 
on the attractive side due to the formation of a spin gap (compare (|6^) and (|59D). For the 
ladder this come from the fact that both in the attractive and repulsive sector, three of the 
modes are always gapped. 

The conductivity above the pinning temperature behaves as 



with again the same exponent than in the d-wave sector ([6^). However, since now Kp^ > 1 
the conductivity now decreases with decreasing T. There will thus be for 1 < Kp^ < 3/2 
a maximum in the conductivity for T ~ Tpin., and the resistivity will go to zero for high 
values of Kp^. This maximum can be seen as a remnant of the superconducting behavior of 
the pure system. For Kp^ > 3/2, the system has infinite conductivity for T — * 0. 



3. CDW sector 

Let us now consider sector IV. In that sector, one has strong fluctuations towards a 
CDW" phase. Such phase is the analog of the CDW^ that existed in the spinless fermion 
problem. We see that the coupling to disorder reduces to (see table |I|) 

j dxCa{x)e'^'+ + H. c. (66) 

As in the spinless fermion case that antisymmetric Charge density Wave only couples to the 
antisymmetric disorder. The RG equation for disorder is 

The antisymmetric disorder is thus relevant for Kp^ < 6. Since the CDW^ phase only exists 
at Kp^ < 1 the CDW^ is always very strongly pinned by disorder. Using ( |67|) we obtain for 
the localization length in that phase 

L^^^-^CDW^ ^{AY^^ (68) 
a \DaaJ 

In the classical limit Kpj^ — > one recovers again the standard resultH for the pinning of a 
classical CDW. 

The conductivity of the CDW^ above the pinning temperature behaves as 

a{T) ~ T^-^ (69) 

showing since Kp^ < 1 a very rapid decrease in the conductivity as T — 0. This behavior is 
a consequence of the very strong pinning of the CDW^. This feature of the antisymmetric 
CDW is similar to the one occuring for the spinless ladder. 
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4- OAF sector 

In the case of the orbital antiferromagnet, the couphng to disorder is made of 2 terms: 
One term comes from Ocdw° the other one from Ocdw^- According to the preceding sec- 
tions, these terms contain respectively cos 0p_ and cos 0o-+ sin 9^^ and (see table ||) there- 
fore have exponentially decaying fluctuations. In order to get nontrivial results, the massive 
modes have to be integrated out as in the preceding sections . This again leads to an action 
of the form (|57| ) and the disorder in the OAF phase is relevant for Kp+ < 3/2. The OAF 
phase is therefore as delocalized as the superconducting SC^ phase, although the pure 
system does not exhibit nay obvious superconducting order parameter. The localization 
length in the disordered OAF is 

Lloc. ^ 

a \Da J 

For Kp+ > 3/2 we have a metallic phase. 

The disorder leads to a conductivity of the form 

a{T) oc t2-2^''+ (71) 

The conductivity in the OAF is therefore identical, as far as the temperature dependence is 
concerned, as the one in the SC^. It will exhibit in the localized phase 1 < Kp^ < 3/2 the 
same maximum in the conductivity for T ~ Tpin.. Once again one sees that the transport 
properties can hardly be guessed from the phase diagram of the pure system. The OAF is 
thus also an excellent candidate for a "superconducting" behavior. 

Using (|36|), it is possible to get some hints on the parameter regime of the extended 
Hubbard model in which the OAF could be achieved. One is in the OAF sector if (^i > 
and Kp+ > 1. In the extended Hubbard language it means 

2V{2 - cos(2A;^a)) < -U < 2V cos{2kpa) (72) 

Let us assumes a local repulsion U > and that one is close to half filling cos(2fci?a) ~ — 1. 
In that case one reaches the OAF for moderate nearest neighbor attraction V < —U /6. Such 
a situation is likely enough to be realized, specially if additional attractive mechanism such 
as phonons are taken into account. 

5. Differences with the spinless ladder 

The spinless ladder and the ladder with spin show some marked physical differences 
. Some of them are due to the fact that interchain hopping has a different impact on 
fermions with spin and spinless fermions. In a system of spinless fermions, energy can be 
gained from hopping only if one site of the rung is occupied and the other one is empty 
due to the Pauli principle. This induces an enhancement of density fluctuations. On the 
other hand, a system with spin can gain energy from interchain hopping by having the two 
sites of the rung occupied by fermions of opposite spins. This leads to spin gap formation 
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(70) 



and a smoothening of density fluctuations. This effect is enhanced in the presence of a 
purely repulsive interaction as it tends to smoothen the density fluctuations in a system 
with spins, whereas it enhances them in a spinless systenS. This has already important 
consequences in the pure case. In particular, the positions of the SC^ and OAF phases 
are different (see figure and figure ^ as the d-wave in the spinless system needs some 
amount of attraction whereas it is achieved from completely repulsive interactions in the 
ladder with spin. In the presence of disorder, the d-wave phase of the spinless system can 
be stabilized by sufficiently attractive interactions, whereas in the system with spin it is 
always unstable(see figures |^ and |^ ), being replaced by an s-wave superconducting phase 
for attractive interactions. Also, in the presence of disorder, the system with spin due to 
the smoothening of the density, shows delocalization compared to the one chain case both 
for the attractive and the repulsive side. On the other hand for the spinless system, the 
reinforcement of the density fluctuations, enhances localization on the repulsive side. The 
attractive side on the other hand is totally delocalized. 

In both case the s-wave phase (occuring for attractive interactions) is very strongly sta- 
bilized by the interchain hopping. This can be understood by a picture of tightly bound 
pairs that behave in both cases as hard core bosons. In that case, the statistics do not influ- 
ence any more qualitatively the transport properties. Similarly both systems tend to form 
charge density waves that are extremely well pinned by disorder (usually much more easily 
than their one chain counterpart). In the case of fermions with spin, this requires some 
mixing of attractive and repulsive interactions so that pair of fermions of opposite spins 
are formed in the chains. These pairs then have hard core bosons interactions so that the 
situation becomes analog to the spinless fermions case. This explains the enhancement of 
pinning for the antisymmetric charge density wave phase. However, in the system with spin 
with purely repulsive interactions there is no CDW^ in contrast with the spinless system. 
Both systems also present an OAF phase that reveals quite stable in the presence of a small 
disorder. For the spinless ladder the OAF is even stable close to the noninteracting point. 
Finally, an interesting similarity between the system of fermions with spin and the system 
of spinless fermions is that pinning on two different CDW phases are possible depending 
on the interactions: either the antisymmetric 2k p CDW or a. 4k p CDW. In these two lo- 
calized phases the behavior of the conductivity at high frequency or high temperature and 
of the localization length at small disorder are very different (the difference appears in the 
exponents) the Akp being much less well pinned than the 2k p- This is to be contrasted to 
the one chain case where only one pinned charge density wave phase is realized. Therefore, 
we may expect to see, for weak disorder, a crossover between two different pinned charge 
density wave phases in the two chain system when varying the strength of the interactions. 
Such a crossover needs a more detailed study. Unfortunately it cannot be tackled by the 
RG since it occurs deep in the localized regime. One interesting question is whether such a 
transition still occurs for strong disorder. 



C. Transport properties 

The ladder with spin shows therefore in presence of disorder transport properties dras- 
tically different from the one one could naively expect form the pure phase diagram. In 
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particular the d-wave phase disappears and does not exhibit any remarkable conductivity. 
Let us look in more details in the transport properties and compare them to what happens 
in a single chainil for the various sectors. 



1. Conductivity 

As was mentioned in section [ill B 2|, the ladder s-wave phase is much more stable to 



disorder than its one chain counterpart (see figures ^ and ^). This effect manifest itself in 
the location of the superconducting-localized transition, and in the localization length. As 
for the spinless case, this effect is entirely controlled by the interactions and going from one 
to two chains affects the power law dependence of the localization length with disorder. It 
is thus much stronger than the increase of localization length occuring for a noninteracting 
system (proportional to the number of channels). In presence of interactions the behavior 
of the localization length cannot be guessed by analogies with the non-interacting system. 
The case of spinless fermions where repulsive interactions make the two chain system more 
localized than the one chain system is an excellent counter-example. 

The resistivity (see (|65D) is also dropping much faster than for one chain for which 
c"ich. (^) ~ T'^^^p. The ladder is thus a much better conductor than a single chain both 
because of the scale of localization and because of the better temperature dependence. In 
addition even in the localized phase the conductivity will increase for all values of Kpj^ for 
which the s-wave phase exists in the pure system, till one reaches the localization temperature 
Tpin.This behavior is qualitatively sketched on figure |^. The SC^ phase shows therefore all 
the "good" characteristics of a "superconducting" phase, and in that respect is much more 
normal than its one-chain counterpart. 

For repulsive interactions a different physical situation occurs. The system is still less 
localized than the one chain counterpart. The transition occurs for a smaller value of 
Kpj^ = 3/2 (versus Kp = 2 for a single chain), and the localization length is larger than for 
one chain (see (|58|)). Contrarily to the single chain where the pinned phase is a random 
antiferromagnet, here the presence of the spin gap forces the localized phase to become a 
pinned ikp CDW. However the SC^ phase is completely wiped out by the disorder, and 
what is more surprising, no trace of this "superconducting" phase can be found in the high 
temperature (T > Tpin) of the conductivity (see (|65D). In particular cr(T) decreases mono- 
tonically even at high temperature in stark contrast with the SC^ phase as shown on figure |^. 
This again illustrate the fact that the transport properties are not linked to the behavior 
of the superconducting order parameter but to the density fiuctuations. For a single chain 
since the density exponent and the superconducting one are related by ^density ~ l/-f^supra 
when superconducting fiuctuation increases density fiuctuations necessarily decrease and 
the system becomes a better conductor. Or course this is also true in presence of a true 
superconducting order in higher dimensional systems. For the single chain the fact that 
superconducting fiuctuations do not necessarily imply better transport also appears one the 
fact that the attractive Hubbard model is more localized than the repulsive one@: when the 
interactions go from repulsive to attractive a spin gap opens and the density fiuctuations are 
suddenly lowered making the system more easy to pin. Similar effect occurs in the d-wave 
phase of the ladder, in a more dramatic way: the d-wave phase do not look superconducting 
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at all since it leaves enough room for enough Akp charge fluctuations. Note that the more 
repulsive the interactions will be the worse the conductivity, in a similar way than for the 
single chain where the phase is a spin density wave. The interchain hopping has thus two 
effects : on the one hand it leads to the appearance of the spin gap that wipes the SDW and 
replaces it by the SC^ wave and on the other hand it freezes the density fluctuations (in 
particular the transverse charge modes). Those gaps suppress 2kF CDW fluctuations, and 
localization happens only through coupling to Akp CDW fluctuations. Since the mechanism 
for localization is the same for all sign of the interactions, the transport properties are only 
weakly dependent of the sign of the interactions. This charge freezing is the dominant effect 
on transport. The two effects are essentially unrelated. 

The most remarkable phase is the OAF which is an illustration of the above. This phase 
has a localization length and a cr(T) as good as a genuine SC^ wave phase ! and yet 
has no genuine superconducting order parameter. In fact the absence of order parameter is 
here also due to the spin gap since for a single chain the corresponding phase is a triplet 
superconducting phase. However the fact that density fluctuations are already very small 
in this phase remains (and is helped by the freezing of transverse charge fluctuations), 
giving the remarkable transport properties of this phase. This remarkable property is not 
an artifact of the potential scattering and persists even if coupling to different form of 
disorder is included. In particular the superconducting-like transport properties of the OAF 
also exist in the presence of a random hopping along the chains and a random interchain 
hopping amplitude (see appendix 0). Note that this phase has analogies of the so-called 
flux phaseii^0, the size of the plaquette is here fixed by the interparticle distance, and of 
course this phase could not be reached for a pure Hubbard model (at the opposite of what 
was claimed in higher dimension). Whether for such a phase, a sort of Meissner effect also 
exists is of course a very interesting question. The connection between the one dimensional 
antiferromagnet and its 2d or 3d counterparts clearly deserves further investigations. In 
particular in two dimensions a phase offering some similarities with the one dimensional 
OAF, has been proposed for the high-Tc superconductorsEi. 



2. Persistent currents 

In a similar way than for the spinless case one can compute the charge stiffness. For the 
ladder with attractive interactions one has for Kp^, Kp < 3/2 

V{L) = V{0) - {—^ f-""" (73) 

Moc.,1 ch. 

whereas for a two chain one, it is 

V{L) = V{0) - {—^ )3-2/^,+ (74^ 

Moc.,1 ch. 

These formulas are valid for a <^ L <^ L\oc.. It is easy to see that they lead to a smaller 
reduction of the conductivity stiffness in the 2 chain case, in agreement with the fact that 
L2 ch. > Lich- For repulsive interactions, it is of the form 
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V{L) = V{0) - (-— ^ f-""^ (75) 

-^loc.,1 ch. 

Whereas for the 2 chain case, it is of the form : 

V{L) = V{0) - {-^ — )3~2i^,+ (7g) 

-^loc.,lch. 

and the 2 chain system has a smaller reduction of conductivity stiffness than the one chain 
system. So up to prefactors the reduction in stiffness in the ladder system with spins is 
identical for repulsive and attractive interactions and the reduction of conductivity stiffness 
also shows no abrupt change as one goes from attractive to repulsive interactions. By con- 
trast, in the one chain case, attractive interactions induce a spin gap and localization arises 
from coupling of a single massless mode to 2kp disorder. This gap closes for repulsive inter- 
actions and localization arises from the coupling of 2 massless modes with the 2kF random 
potential. This causes the abrupt change in transport properties and charge stiffnesJilS 
when one goes from attractive to repulsive interactions. This is related to the fact that the 
localization lengths for attractive and repulsive interactions have the same dependence on 
disorder, in marked contrast both with the spinless problem and the single chain with spins. 
The effect of increase of persistent current by repulsive interactions occuring in the single 
chainEj is thus either absent or strongly reduced (not an exponent effect any more) in the 
ladder. It would of course be interesting to investigate in ladder with more than two legs 
to see if this effect reappears and check for possible difference of behavior between odd and 
even legs ladders. 

IV. COUPLED LADDERS 

A. Mean field treatment 

In the preceding sections, we have been considering isolated bichains. To describe realistic 
compounds, such as SrCuO, and have a finite temperature phase transition, interchain 
coupling should be taken into account. A realistic coupling is of course single particle 
hopping between the ladders. However in ladders, due to the existence of single particle gaps 
(spin and antisymmetric charge mode) for the ladder, single particle hopping is irrelevant, 
provided that the inter-ladder hopping is much smaller than the gaps of the system. One 
has therefore to consider onlv the particle-hole (or particle-particle) coupling generated 
by the single particle hoppingp. Such couplings can lead to an ordered phase at a finite 
temperature. As is very reasonable on physical grounds such interchain couplings stabilize 
the dominant one-dimensional instability. We focus here on the existence of a stable d-wave 
superconducting phase. This allows us to keep only the particle-particle (or Josephson) 
coupling between the ladders. The Hamiltonian for the coupled ladders system is: 

J 



H — ^2 -^disordered 2 chain system.n + 2" / ^•'^^^SC,n{^)0 SC,n+l{x) + Olc:n+l{^)OsC,n{x)) 



(77) 
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Where Osc,n is the operator for (d-wave or s-wave) superconductivity for the n— th ladder 
and J is the strength of the Josephson couphng. On can simphfy further the Hamihonian 
(fr?!) by keeping only massless modes in the ladder. Doing so we assume that the spin gap and 
the interchain gap of the two chain system are much larger than the disorder and much larger 
than the Josephson coupling. However, we make no assumption on the relative magnitude of 
the Josephson coupling and the strength of the random potential. The resulting Hamiltonian 
is, both for the case where the dominant instability is s-wave or d-wave superconductivity 



2n 



dx 

+ 1 — 



ecfr.,n(a:)e*'^''+-+H. c. 



+ J J dx COs{9p+^n- Op+^n+l)) 



(78) 



To solve (^) we treat the Josephson coupling in mean field assuming the existence a fi- 
nite superconducting order parameters (cos(6'p+)). By making the replacement cos(6'p+^„ — 
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p+,n+l) 



(COS 



(^P+))cos(^, 



the Hamiltonian (ffSf) becomes the one of an isolated ladder 



system in an external field, the value of which is determined by a self-consistency condition. 
The Hamiltonian is then 



H 



MF 



dx 
W 
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Up+Kp+{TiT[p+f + 



+ 



dx r 



U{x)e'^'^'^ +Y{. c. 



(27ra)^ 



dx cosf 



(79) 



with the self-consistency condition W = J{cos{9p+)). 
The equation determining is 



1 1 /• /•'Sc 



(80) 



with = Y' -^0 is Hmf for W = 0. To solve (|80|), one has to compute the finite temper- 
ature superconducting response function of a ladder in the presence of disorder. There are 
presently no methods to do this exactly, but one can get an accurate solution for Tc by mak- 
ing some simplifying approximations^. First, one notices that a finite temperature induces 
a cutoff length 1{T) = Up^/T beyond which all correlation functions decay exponentially to 
zero. We make thus the approximation that beyond 1{T) all correlation functions are truly 
zero and below 1{T) they are equal to the T = correlation functions. This allows us to use 



the RG equations introduced in section |T|. If we denote by x the superconducting response 
function, when we change the running cutoff a{l) a{l)e'^^ we have x ~^ 2k\i) ^- 
Thus to compute correlation functions at lengthscale R it is sufficient to integrate the RG 
equation from the cutoff up to R and follow the renormalization of the response function. 
Making use of these two approximations, the equation giving Tc simplifies into : 



1 
J 



RdR 

^exp( 



HR/») dl 



27ra2"--v 2K{1)' 
The values of K{1) are obtained by numerically solving the RG equations: 



(81) 
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— = -Dil)Kiir (82) 
^ = (3 - 2K{l))D{l) (83) 

the values of forK = 0.5, 1.2 and J = 0.1 as a function of D are shown on figures ^ and 
H respectively. We note that for = 1.2 we have an s-wave superconducting phase and 
for K = 0.5, a d-wave phase. This can be expected since the interchain coupling stabilizes 
the dominant one dimensional fluctuation fsee figure We see that (see fig. |]) as in the 
case of the single chain mean field theoryO of superconductivity we have an initial linear 
decrease of the critical temperature with disorder strength. This is to be contrasted with the 
standard mean field theory of the s-wave superconductor in three dimensions being based on 
a diffusion approximation that does not include Anderson localization effects and gives Tc 
independent of the disorder. This is the well known Anderson theorem. The linear decrease 
of Tf. with the strength of disorder for s-wave superconductivity in our chain mean field is due 
to localization effects. This peculiar situation is due to the absence of a diffusive regime in 
one dimensional disordered systems, which implies that their response functions are always 
affected by localization effects. 

For a d-wave superconductor one expects in mean-field theory a linear decrease of Tc as a 
function of D (see e.g. Ref. |7^). For the ladder system however the decrease of Tc is mainly 
due to the localization effects, similarly to the s-wave superconductor case. Although it 
indeed starts linearly for small disorder (see figure 0), localization effects manifest themselves 
by the the sudden drop to Tc = at a critical disorder strength (see fig. ||,0). 

For identical Josephson coupling between the bichains, the critical disorder strength is 
smaller for the d-wave superconductor than for the s-wave one. 



B. Simplified Treatment: 

Although the mean field theory allows an accurate description of the effects of disorder 
on Tc the critical value of disorder above which superconductivity is destroyed can also be 
obtained by a very simple physical argument. Let Tj:^^^'^\j) be the temperature at which the 
superconducting transition would occur in the array of ladders if there were no impurities. 
Just above T^p^^'^)^ the thermal length is ^^urc) and beyond that length all phase coherence 
is lost. Clearly, if the thermal length is smaller than the localization length in a single 
chain containing impurities ^loc, phase coherence is lost before coherent backscattering can 
build Anderson localization. The system will escape localization due to the building of the 
(mean-field) superconductivity. Thus, if 

""'^ <6oc (84) 



rpipuve) 



Anderson localization will not suppress the superconducting transition. Equation (|^) gives 
a simplified criterion for the stability of superconductivity. For fixed Josephson coupling 

J, (|80|) leads to T^p^'^'^) ~ J^''p+-\ Thus, the higher Kp+ the higher tJp"^^). From the 
preceding section, the localization length both in the s-wave and d-wave superconducting 
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phase is ^loc. ~ (jy)'^'^'^''^ ■ Increasing Kp+ also reduces ^loc- Thus, the two effects reinforce 
each other, and make the SC^ phase, that exists for Kp+ > 1 more stable against Anderson 
localization than the SC^ phase that exists only at -ft'p+ < 1. 

V. EXPERIMENTAL CONSEQUENCES 

The theoretical results obtained in the preceding sections have important consequences 
for experimental systems that are believed to be well approximated by coupled chains sys- 
tems, namely the doped ladder systems which present a superconducting transition and the 
2 band quantum wire. In the former case, one would like to know if the superconducting 
transition is related to the divergence of superconducting fluctuations in the strictly one di- 
mensional system that results at the mean field level in a finite T superconducting transition 
or if the physics of the transition is a two or three dimensional one. We believe that the 
resilience of superconductivity to disorder is a stringent test of effective dimensionality. In 
the case of quantum wires, we discuss the experimental consequences of our results for the 
conductivity and charge stiffness in the interacting system. Measurements of the conduc- 
tance would allow to check th above theories for the ladders and provide a measurement of 
the Luttinger liquid parameter in the charge sector, providing some insight on the strength 
of interactions in these systems. 

A. Superconductivity of doped ladder systems 

Our study has various experimental consequences for the observation of superconduc- 
tivity in Ladder systems. First, if the superconductivity is to come from purely repulsive 
interactions (i.e. to be of the d-wave type), it should be extremely sensitive to disorder as 
we showed in section [IV A| . In fact, any randomness would induce a conductivity that never 
increases as temperature decreases (see figure ||), so that superconductivity would be impos- 
sible to probe except in extremely pure samples. Such sensitivity with respect to disorder 
is certainly consistent with the difficulty in observing any type of superconductivity in the 
ladder systems Srn_iCuri4-i02n^^- However superconductivity seems indeed to be observed 
in Sro.4Cai3.6Cu2404i.36lIj under pressure (~ 3GPa). Whether such superconductivity is of 
the d-wave type is of course still open. Various experimental facts, however seem to indicate 
that if it is the case, it is unlikely that such superconducting phase could be described by 
weakly coupled ladder systems. Indeed one could use the criterion ( p4|) to estimate the lo- 
calization length. Taking a reasonable value of 10^ms~^ for the Fermi velocity, one obtains 
from the observed ~ lOK, a minimal localization length of ^ ~ lOOOOA. Using (p8D, this 
leads to extremely long mean- free paths {I = for K = 1/2) when one is in the d-wave 
sector. So unless the chains are extremely pure, a fact not likely to be true in such doped 
materials, one expects based on one dimensional physics alone that the superconductivity 
should be totally suppressed. If the presence of superconductivity is due to an extremely 
pure system (which is doubtful) then, introducing more disorder in the system (for instance 
by irradiation) should induce a dramatic decrease of the critical temperature. 
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Besides the extreme sensitivity of Tc to disorder other arguments are again a simple 
stabihzation of one-dimensional physics in the experimental compound: even if one could be 
below the critical disorder strength determined by ( p^ and Figure 0, the physics above T^. 
should be dominated by the one-dimensional (ladder) effects. In this regime the resistivity 
goes up with decreasing temperature as described in section |111 C| . The observed resistivity 
showing a monotonic decrease of the resistivity (roughly with a or T law) is again 
incompatible with the one-dimensional description. If one is in the purely repulsive sector, 
the most likely explanation of the main experimental features is that under pressure the 
interchain hopping between the ladders become strong enough so that the system does not 
retain its one- dimensional feature, but is more accurately described by two-dimensional 
physics. Such an interpretation is also compatible with the fact that the system at ambient 
pressure is insulating. In that case, the coupled bichains treatment becomes extremely 
questionable, and it is probably better to start from a two-dimensional description, for 
which disorder effects are probably weaker, but for which the nature of the superconducting 
phase has yet to be completely elucidated. 

Another interesting, but probably more farfetched, possibility could be that the system 
is in fact in the orbital antiferromagnetic sector. In such a sector the effects of disorder 
are much more reduced, and even very large localization lengths can lead to reasonable 
mean-free paths (^/a is at worst {l/aY, and diverges for i^" = 3/2, see (|T0|)). The resistivity 
decreases with temperature according to (|7TD . Here the difficulty lies more in getting the in- 
teractions corresponding to this phase, since one needs local repulsion and a sizeable nearest 
neighbor attraction. In any case, careful measurements of the temperature dependence of 
the resistivity above Tc could help to decide if such OAF effects are present. Of course here 
again, one cannot exclude that the physics is two-dimensional to start with, but at least 
now the one- dimensional starting point is more consistent with the dominant experimental 
features. 



B. Application to quantum wires 

Progress in nanostructure technologies have allowed for measurements of the transport 
properties of low dimensional electronic systems. In particular, in recent experiments on 
quantum wires^aEl, the conductance of a quasi one dimensional electron gas has been 
measured at very low temperatures. For the pure system, or extremely weak disorder 
one finds quantized values of the conductance HTI in good agreement with the theoretical 



predictionsO'El at fractions of ^ as a function of width of the quantum wire (i.e. of the 
number of subbands at the Fermi level). The relation between the number of channels and 
resistance has been verified^. Impurities on the other hand induce backward scattering 
that is known to cause Anderson Localization in a sufficiently long system. In small enough 
system, it leads to reduction of conductance as the length of the system is increased or the 

2 

temperature is lowered. Deviations of conductance from ^ as a function of temperature 
have indeed been obtained in experiments^ as well as deviations as a function of of the 
length of the wire0 and can be related to the Luttinger liquid exponent. The correction to 
the conductance due to impurities&ii is of the form 
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G{T) = j- gT-^ (85) 

where z/ = 1 — Kp is the conductivity exponentil. The derivation is similar to the derivation 
of the T dependence of conductivity in appendix For finite size systems T can be replaced 
by the lower cutoff Vp/L in (|85D. This formula only holds at high enough temperatures or 
for systems of length L shorter than the localization length for which the corrections term 
is small. 

If two channels are present in the wire, the system becomes then equivalent to a lad- 
der system. Two band present at the Fermi level are the equivalent of the bonding and 
antibonding bands of the ladder system. On then expect that the whole physics derived 
in section |111 C| should apply to these wires. In particular, since one expects reasonably 



repulsive interactions one should be in the SC^ of the CDW'^^'' phase. Going from a single 
chain to the ladder should have observable consequences on the transport properties . First 
since the localization length increases drastically in the ladder system one would expect the 
conductance corrections due to disorder, to be much weaker for two channels. This of course 
assumes that the typical interactions do not vary too much when going from one channel to 
two channel, a fact which is not certain. Secondly by doing the expansion of the corrections 
to conductance for the ladder and using the conductivity exponent ( |1II C] ), one would obtain 
u = 2 — 2Kp^. A fit of the temperature dependence of the conductance in§, could allow 
to extract the Luttinger liquid exponents for the ladder system, as well as check the above 
predictions^. 



VI. CONCLUSION 

In this paper we have examined the effects of disorder on a 2-legs ladder system, using 
RG techniques. We have computed the effects of disorder on the phase diagram as well 
as the localization length. Disorder has drastic effects on the phase diagram. For spinless 
fermions, it leads to an extremely strong localization of the charge density wave phase 
that exists for repulsive interactions. Such localization is even stronger than for a single 
chain. On the other hand for the ladder system there is a remarkable stability of the s-wave 
superconducting phase (for attractive interactions), compared to the single chain case. The 
insulator-superconductor transition occurs in the vicinity of the non interacting point for a 
pure t-V model whereas in the one chain system it occurs for strongly attractive interactions. 

For fermions with spin, the repulsive part of the phase diagram is also strongly localized 
by disorder. In particular the novel d-wave superconducting phase found for ladder systems 
is completely suppressed by an arbitrarily small amount of disorder. We emphasize that this 
is not only a pair breaking effect but a much stronger Anderson localization effect. On the 
other hand the s-wave superconducting phase occuring for attractive interactions is again 
much more stable to disorder than its one chain counterpart. 

Besides obtaining the phase diagram, we have also investigated the transport properties 
of the t-V and Hubbard two chain systems. The RG enabled us to compute the localization 
length and the charge stiffness as a function of disorder (see tables fT|,fVp and the tempera- 
ture and frequency dependence of the conductivity. Various remarkable fact emerged. First 
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the behaviors of the spinless ladder and the ladder with spins are very different. In par- 
ticular the spinless ladder shows the same tendency than the single spinless chain, namely 
that attractive interactions decrease localization whereas repulsive interactions enhance it. 
In the two chain case, that effect is even stronger. For attractive interactions there is no 
localization, whereas for repulsive ones the system is much more localized than its one chain 
counterpart (the exponent in the dependence of localization length with disorder is changed). 
On the other hand for the Hubbard ladder there is no such effect: up to a prefactor, the 
localization length with attractive interaction is the same as for repulsive ones. As a conse- 
quence, the corrections to conductivity stiffness are the same for attractive and for repulsive 
interactions. 

The temperature dependence of the conductivity follows a power law of the form o"(T) ~ 
rp2-2Kp+^ for temperatures above the localization scale Tioc, where Kp. For the repulsive side, 
where for the pure system one would have the d-wave superconducting phase, Kp+ < 1, and 
thus the conductivity decreases as a function of the temperature even well above Tioc. The 
transport thus shows no remnant of the superconducting behavior one could have naively 
expected when looking at the pure system. This remarkable fact illustrates that transport is 
in fact controlled by the density fluctuations of the system and not by the existence of not of 
slowly decreasing superconducting correlation functions. Ladder system provides evidence 
of a phase that is genuinely a d-wave superconductor as far as phase diagram is concerned, 
but would from the transport point of view be closer to an insulator. Of course such an 
interesting behavior would clearly deserve more studies. In particular it would be interesting 
to know how the correlation between the density fluctuations and the superconducting one 
evolves as the number of chain is increased, and how the crossover to the three dimensional 
situation occurs. Such a study goes of course far beyond the goals of the present paper. 

We have applied our results to two types of experimental systems. First our results should 
be relevant for quantum wires with two channels. Here the prediction for the exponent in 
the conductivity can be directly checked by the measuring temperature dependence of the 
conductance of the system. Note that the conductivity/conductance exponent 2 — 2Kp+ for 
the ladder systems is different from the one for a single channel (or a single chain) 1 — K. 
Due to the increase of the localization length when going from one channel to two channels 
one would also expect overall smaller corrections to the conductance for a given strength 
of the disorder, and roughly constant interactions. Investigation of systems with more than 
two chains would be useful in order to get a better understanding of the role of internal 
symmetries and gaps on the transport properties of quasi one dimensional systems. This 
is of course also useful in connection with experiments on quantum wires. In particular, 
we expect that the behavior of systems with an even number of legs is dominated by gap 
formation whereas the behavior of systems with an odd number of legs should be closer to 
the one of a one chain system. 

The other experimental system on which our results could be applied are of course the 
coupled bichains where superconductivity has recently been obtained!^ under pressure. To 
compare with this system of coupled ladders, it was necessary to treat the coupling between 
different ladders, which we did using a mean-field approximation. The results indicate that 
even in the presence of coupling between ladders the d-wave phase is still much too sensitive 
to disorder, to be the one experimentally observed. In addition the observed temperature 
dependence of the conductivity would be incompatible with the one computed here should 
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these systems be dominated by one-dimensional (ladder) physics. These observations and 
the fact that the conductivity occurs under very large pressure tends to indicate that the 
mechanism for superconductivity in these systems is very likely to be of a two or three di- 
mensional nature and not just the mere stabilization of the ladder superconducting phase. 
On the other hand the system without pressure has a resistivity that could be more com- 
patible with the localization effects described here. Of course one interesting question would 
be whether one can get a one (ladder) to higher dimensional crossover as the pressure is 
applied. This of course could only be decided by a more quantitative comparison with exper- 
iments as well as further theoretical and experimental work. Adding additional impurities, 
for example by irradiation, could allow to distinguish if the system is in a one-dimensional 
regime, since one expect much more drastic localization effects in that case. 

Finally a ladder system with spins exhibits an extremely interesting orbital antiferro- 
magnetic phase. Although such a phase cannot occur in a pure Hubbard system it can in 
principle be stabilized if some nearest neighbor attraction is added. Although such a phase 
has no superconducting order parameter, it has perfect conductivity in the presence of a 
random potential. Moreover that perfect conductivity is also robust in the presence of ran- 
dom hopping both along the chains and perpendicular to the chains. As far as transport is 
concerned this phase is therefore a one dimensional "superconductor" . Nevertheless, it has 
only subdominant (in the spinless case) or exponentially decaying (in the case of fermions 
with spin) superconducting correlations, again an illustration that looking at the supercon- 
ducting fluctuations is not a good criterion to determine the transport properties. Due to 
the peculiar nature of this phase it would be interesting to check whether it survives in lad- 
der systems with more than two legs. More generally it also deserves further investigation 
in dimensions higher than one, both in relation of flux phases of two dimensional systems 
and other orbital phases proposed for the normal state of cuprate superconductors!^. 

The study of the disorder effects could also be extended in various directions. In partic- 
ular a more detailed description of the physics inside the localized phase would be suitable. 
However such a description is beyond the reach of the simple RG calculation. Going to 
strong but diluted disorder is also a challenging problem. In particular understanding the 
crossover from the results of our paper to the limit where disorder suppresses gaps altogether 
in the system remains yet to be done. 
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APPENDIX A: BOSONIZATION TECHNIQUE 

In this section, we will give a short review of the bosonization technique in order to fix the 
notations. We give the expressions for a single chain of spinless fermions. For more species 
of fermions, one can bosonize each specie individually, and the corresponding expressions 
are given in the text. 
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1. representation of fermion operators in terms of boson ones 



Non-interacting one dimensional spinless fermions on a lattice are described by the kinetic 
energy 

N 

H = -tJ2{cl^^Cn + cicn+i) = ^ e(A;)4cfc (Al) 

n=l k 

where e{k) = — 2tcosfc and c„ = -^J2k^^'^^"'^k- 

To obtain the asymptotic (low energy, long wavelength) properties of the system one 
can linearize the spectrum near the Fermi "surface" {^kp) and take the continuum limit by 
introducing ijj{x) = ^ with a the lattice spacing and x = na. With our definition the ip^s 
have the commutation relations of continuum fermion operators. We define the R (resp. L) 
(right and left movers) fermions as fermions with momentum close to +kF (resp. —kp) as 

M^) = E e'''c,^+k (A2) 



and similarly for ipiix) with kp ^ —kp- A is a cut-off needed not to double count fermion 
states, and imposed by the linearization of the dispersion relation. All asymptotic properties 
can be expressed in term of In particular the full fermion operator becomes ipi^) = 

e^''''^ipji{x) + e~^^''^ipL{x). The Hamiltonian (|A1|) becomes 

H = -ivp{ilj'^^d^iljR - ^ida^^l^L) (A3) 

with vp = 2tasin{kpa). 

Due to the separation into two branch of fermions and the linearization of the spectrum, 
the Fourier components of the fermion density operators 

Pb.,l (q) = E '^R,L,k+q^R,k ( A4) 

k 

have boson commutation relationdli'S 

Ztx 

[pL{q),pL{-q')] = ^(l^m' 

[pL{q).PR{-q')] = Q (A5) 

This allows to rewrite f|A3) as 



H = iTVp J dx Pr{xY + Pl{xY (A6) 

with Ps(x) = ip\{x)'ips{x) for s = L,R. Instead of using the density operators themselves it 
is more convenient to introduce 
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n(x) = PR- PL 

—dA = {pR + Pl) (A7) 

TT 

Physically, 11 is a momentum density while dxcf) is proportional to the deviation of the 
fermion density from its average value. The commutation relations for the p's imply that 
n(x')] = i5{x — x'). Also, the Hamiltonian rewritten in terms of 11 and is : 



H = I dx— ^ 

271 



[711 



+ {dx4>f (A8) 



which is just the continuum limit of the Hamiltonian of a ID harmonic chain. Note that the 
following procedure could have been applied to a more complicated lattice Hamiltonian than 
(|AT|). All that is needed is that the Fermi surface reduces to two points. The effectiveness 
of bosonization stems from the fact that it is possible to express the fermions operators 
in terms of H(x) and If one introduces 6{x) = 7c J^^Il{x')dx' one has the following 

relations: 



V'l(x) = — ^e*[^W+<^(-)l[/^ (A9) 
y/27ra 

a being a cutoff the presence of which is imposed by the cutoff needed in the linearization 
of the dispersion relations. 

the Ur and Ul are anticommuting operators introduced by Haldane that annihilate one 
fermion at the Fermi level. These operators also anticommute with their hermitian conju- 
gates. It can be verified explicitly that those relations reproduce correctly the commutators 
of fermion operators. These U operators give in general corrections vanishing in the thermo- 
dynamic limit and can be safely dropped. On the other hand, if there are different species 
of fermions (such as up and down spin fermions or band degeneracies), one must bosonize 
separately each fermion specie using the formulas for spinless fermions. It is needed to in- 
troduce UL,n,UR^n operators and their complex conjugates (n indexing the internal degrees 
of freedom such as spin) to enforce proper fermions anticommutation relations. In order to 
make that bookkeeping less tedious^, one can introduce r] operators such that : 

Vi = Va (AlO) 

where a = {L, n) or a = {R, n) these operators can therefore introduce minus signs in the 
various bosonized expressions. 



2. handling the interactions with bosonization 

Let us consider spinless fermions. Interactions can then be handled straightforwardly: If 
one adds a density density coupling of the form / dxU p{xY The density can be decomposed 
in a slowly varying part Pr{x) + Pl{x) and a 2kF part e^'^^''^'ilj\{x)tljR{x) -\- H. c. 
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In the Hamiltonian, one retains only the slowly varying terms (the other term give a zero 
value when integrated over x). The 2kF always disappear, while the Akp can persist in a half 
filled lattice system@. As a consequence, at a non commensurate filling, the Hamiltonian 
reduces to: 



H= — 

271 



K 



(All) 



with uK = f as a consequence of galilean invariance. if one makes the rescaling (p (p / y K 
and n — > Y{\/K one has the same Hamiltonian as in (X8) with the correct commutation 
relation for and H. If one computes the physical correlation functions at OK such as 
the 2kF part of the fermion Green's function G{x — x',t — t') = —z{TipR{x,t)'ilj\^{x',t')) it 
is easily seen that K controls their power law decay while u controls the propagation of 
excitations, u and K are also related to physical quantities such as the charge stiffness!^ 
and the compressibility. More specifically, defining the compressibility by x = —^^^^ 



P = -(If )t and taking T 



OK, we have : 



X 



ttK 

ukp 



(A12) 



The charge stiffness is defined by T> 
^From that definition, one obtains : 



^{^-^^)^=Q) ^ being a fiux threading the system. 



V = uK 



(A13) 



The case of fermions with internal degrees of freedom is usually more complicated, be- 
cause some of the interactions cannot be reduced to {dxCpY terms, the most well known 
example being the backscattering of two fermions with opposite spinsi'i. Usually, one finds 
sine- Gordon Hamiltonians of the form: 

dx 



= ' 2i 



u 



uK{'KnY + -{d, 



+ A J dx cos(/?0) 



(A14) 



These Hamiltonians can be studied using renormalization group (RG) techniquesBl^. The 



A has scaling dimension 



fiow equations for K and A are of the Kosterlitz-Thouless form0 ! 
2 — f3'^K/A. therefore a small A is relevant for K < 8//?^. /^From the RG equation for A 
one sees that there are two regimes: one small K or large enough A regime, where A is 
relevant and a large small enough A regime where A is irrelevant. When A is irrelevant, 
the correlation functions keep their power law character up to logarithmic correctiong^a. on 
the other hand, if A is relevant, will acquire a non- zero expectation value that minimizes 
the ground state energy and a gap will be generated. It can then be shown! that there 



(/(0)) ~ fm) and that {T^e 



exp^ 



-) where ^ is a correlation 



length. These results are used extensively in the paper. 



APPENDIX B: MEMORY FUNCTION CALCULATION OF AC AND DC 

CONDUCTIVITY 

For the sake of clarity, we will explain the technique on the example of one chain of 
spinless fermions (technically this is the simplest case) , and then explain how the calculation 
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can be extended to more complicated cases. First let us describe the memory function 
approximationii. The conductivity is given by linear response theory as: 

aiu) = _,^W^lM (Bl) 

LJ 

where x{^) is the current- current response functioniiS. The memory function M(a;) is 
defined by: 

= (B2) 

This gives the exact formula!^ : 

M(.) = (B3) 

an expansionii at high frequency and small impurity concentration gives : 

_ («F;F».-<<f;F»..o)/c. 
-X(0) 

Where -C; ^ is a retarded correlator evaluated for the pure system and F = [J, H] J being 
the total current. To use that formalism in the framework of bosonization, we first need 
an expression for the currentS. This can be obtained from the definition of the fermion 
density p{x) = Pq — ^ and the current conservation equation : dtp + d^j = 0. One obtains 
j{x) = —. Using the Heisenberg equation of motion for and noting that the total current 
J = J dxj{x) on finds J = uK ^'n.{x)dx. The coupling to disorder being: 

i/imp = / dx^e'"^^^^ + H. c. (B5) 

We get F oc / |Me*2<^(-) - me-*2<^(-). 
This gives : 

(T,F(r)F(0)) oc / dxD6{x){ „ ] oc r-^^ (B6) 

Therefore, < F; F >^oc J dre'^^'^T-'^^ oc cu^^^^ This gives M{uj) oc cu^^f-a ^nd for 
K > 3/2 (j{uj) oc uj'^~'^^ . The formula we have obtained is valid only at high frequency. 
We can get from it a high temperature formula by using the dimensional equivalence of 
temperature and frequency ( e. g. hu ~ ksT). To generalize the calculation to a more 
complicated case, we must first note that in the formula for the current, will be replaced 
by (f)p in the case of two chains of spinless fermions and 0p+ in the case of 2 chains of fermions 
with spin. The coupling to disorder being some / dx^{x)e^^^ + H. c. , n depending on the 
problem at hand, we see that in the general case we will just have to make the replacement 

2 

2K — > ^^K in the formulas giving (T(to'), ct(T). 
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APPENDIX C: EFFECTIVE RANDOM POTENTIAL IN THE PRESENCE OF 

GAPS 



In that section, we will give a derivation of the RG equation for Da at > 
0. We start with the method oiS : we compute perturbatively the correlation func- 
tion : ^Xre«v^</'p(3::i:n)g-«v^</>p(2^2,T2)^_ jj^ second order in the random potential, since 

(r^sin(\/20||)(x,r)sin(V20|j)(O,O)) ~ e-*^/', there is no singular contribution. Therefore, 
we must go to fourth order. We will drop the combinatorics since we are only interested in 
the renormalization of D. The fourth order term is of the following form : 

2 f dXidTidX2dT2dx^dT2.dXAdTAr^, . 
L)„ / 7 77 [o[Xi - XA)b[X2 - Xs) + d[xi - X2)o[xz - Xi)\ 

jrp ^l\/2[<f>p(x,T) + <f>p{xi,T-i_) + <l)p{x'J,,Ty,)-<l>p{x2,T2)-<t>p{x4,,Ti)-(f>p(Qfl)]^ 

{Tr sin(V20||)(xi, Ti) sin(V20||)(x2, ts) sin(y20||)(x3, rg) sin(V20||)(x4, u)) (CI) 

The 0|| will be exponentially small except when | ri — |-C / and | r2 — |<^ / or 
I '"i — ''"2 l-C / and I rs — r4 |-C / (the other cases are equivalent to these two ones up to a 
relabeling of dummy integration variables). It is easily seen that the second case is in fact 
trivial. Therefore, the only interesting contribution comes from the first term. This term 
reduces to the simple form : 

DlfC I dXidndX2dT25{Xi - X2)(T,e*^[^''("'")+2'^''("i'"i)-2*''("2,r2)-^,(0,0)]^ ((^2) 

Where C is a constant that depends on the regularization of the continuum model. It can 
be seen that the term that we obtain can be generated by the following effective coupling : 

//effective = / dxU{x)e''^'^^ + H. c. (C3) 



with ^eff.(a^)^eff. (^') = D6{x — x') and D oc CD^. It is clear that for couplings of the 
form e*'^'' cos(^|| ) the same argumentation will be equally valid. Note that using SCHA 
approximation gives different results; This is due to the fact that normal ordering in SCHA 
is done without taking the presence of the gaps into account. Therefore standard scaling, 
irrespective of the presence of the gaps always holds when one uses SCHA. 

APPENDIX D: THE ORBITAL ANTIFERROMAGNET IN THE PRESENCE OF 
RANDOM INTRACHAIN HOPPING AND RANDOM INTERCHAIN HOPPING 

We consider the following two types of random hopping: a random hopping along the 
chains: 

-f^intrachain = X! [^^i (4+l,(T,lCi,<T,l + c1,ct,1'^«+1-'^.i) + ^^1 i'^l+l,a,2Ci,cr,2 + 4,^,2^1+1,^,2) (d) 
i,cr 

a random interchain hopping amplitude: 
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-f^interchain — X] '^'^-L,* ('^j,l'^«,2 + Cj_2Ci,lj (D2) 

Where 5t± is rea/. Bosonization of (Pl|) leads to an expression identical to the one that 
obtains by bosonizing a random on-site potential. It is then evident that the transport 
properties of the Orbital Antiferromagnet are the same in the presence of a random potential 
or random hopping along the chains. Bosonization of equation ( |D2| ) gives the following 
expression: 

-^'t'^" {x)e~'-'^''+ [i sin (/)p+ cos cos + cos (j)p- sin sin 0^+] + H. c. 

(D3) 

It is not difficult to see that such term has exponentially decaying correlations since 9p_ 
develops a gap. Integration of the massive p— mode leads to a coupling that is identical 
to the coupling to a random potential. Therefore a random amplitude of the hopping term 
also does not affect the transport properties of the OAF more severely than a random po- 
tential and thus the "superconducting" transport properties of the OAF are not an artifact 
of restricting to random potentials. All physically admissible random perturbations of the 2 
chain system lead to the same limit for localization delocalization (-ft'p+ = 3/2) the same be- 
havior for conductivity as a function of frequency and temperature and the same dependence 
of localization length as a function of disorder. 

On the other hand if the random hopping term has a random phase, there is a direct 
coupling to the OAF order parameter and then the OAF phase is suppressed. Such terms 
are allowed for instance in a tight binding picture only if the phases on the atoms of the 2 
chain system cannot be made real. This could be achieved with a random magnetic flux in 
each plaquette of the 2 chain system. 
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TABLES 



TABLE I. The 4 sectors of the pure 2 chain spinless fermions model, as a function of Kp and 
gf. The average value of the massive field (^y) are indicated together with the phase with the most 
divergent susceptibility 



I 


II 


III 


IV 


9f + 
Kp < 1 


+ 

> 1 


> 1 


< 1 


(^ll) = ^ 


TT 








phase OAF 






CDW 


TABLE II. The 4 sectors of the pure 2 chain Hubbard model, as a function of Kp and gi. The 
average value of the field developing a gap are indicated together with the phase with the most 
divergent susceptibility 


I 


II 


III 


IV 


91 + 
Kp+ <1 


+ 

> 1 


> 1 


< 1 


(Op-) 

f 

phase SC^ 




TT 

{0a-] = 

OAF 





=0 





(Oa-) = f 

CDW 


TABLE III. The conductivities and localization lengths in the spinless fermions case. The 
phases are the ones of the non disordered system that are turned into localized ones upon intro- 
duction of a small disorder. 


phase 






a(T) 








rp2-4Kp+ 


PCDW 


(At)"-"'* 




rj'2-AKp+ 


TABLE IV. The conductivities and localization lengths in the fermions with spin case. The 
phases are the ones of the non disordered system that are turned into localized ones upon intro- 
duction of a small disorder. 


phase 






a{T) 


OAF,SC'^,SC' 






rp2-2Kp+ 


CDW 


ik)'-''* 
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FIGURES 



OAF 



SC 



K 



CDW 



SC 



FIG. 1. The phase diagram of a generic spinless ladder in terms oi gf and Kp. Kp > 1 means 
attraction in the symmetric charge sector and Kp < 1 repulsion. The line depicts the phase spanned 
by the pure t-V ladder, leading to a CDW^ phase for F > and a superconducting SC^ phase 
for F < 0. 



38 



FIG. 2. The phase diagram of the disordered 2 chain t-V model in terms of gf and Kp. For a 
single chain the system is localized for Kp < 3/2. Ladder effect therefore delocalize for attractive 
interactions and enhance localization for repulsive ones. 
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4k 

CDW ^ 


SC^ 







OAF 


1/2 




1 


K 

P + 


CDW" 






sc 



FIG. 3. The phase diagram of the pure 2 chain Hubbard model in terms of gi and Kp^. Kp^ > 1 
and gi < corresponds to purely attractive interactions. Kp^ < 1 and gi > 1 corresponds to purely 
attractive interactions. For a Hubbard model, this leads to a SC^ phase for i7 > and a SC^ 
phase for [/ < 0. 
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/ / / /^CDV\ 


/ / / 

A\£ / / 
/4K / / 

i y / 


OAF 




A / / 

/ / 4k / 

PCMjW / 

■ WW VV / 


3/2 Kp ^ 

SC^ 



FIG. 5. The phase diagram of the disordered 2 chain Hubbard model in terms of g\ and Kp^. 
The Sd^ phase is completely eaten by the PCDW^^^ phase whereas the OAF and the SC^ persist 
if there is enough attraction. Delocalization occurs for K > 3/2 i.e. for less attractive interactions 
than in the one chain case. 
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a(T) A 




loc 



FIG. 6. Behavior of the conductivities of the s-wave (dotted hne) and d-wavc (sohd hne) su- 
perconductor as a function of temperature. For T 3> Tioc., (^{T) oc r2-2i4"p+ the d-wave, there 
is no maximum in the conductivity and therefore no remnant of superconductivity in the locahzed 
phase. 
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0.6 - 
0.4 

0.2 - 



0.05 0.1 0.15 0.2 0.25 0.3 

D 

FIG. 7. Tc as a function of disorder for the d-wave phase {Kp = 0.5). Tc drops quickly to 
D ~ 0.2 



44 



Tc 



1 •*•■ 




0.4 - 



0.2 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 8. Tc as a function of disorder for the s-wave phase {Kp = 1.2) . Tc drops to zero for 
D ~ 0.9. Note the initial hnear decay of Tc that shows that Anderson Theorem does not hold in 
coupled chain system due to strong localization effects and absence of a diffusive regime. 
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